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Abstract. In [8, §8], the first author outlined an algorithm for calculating a parametriza- 
tion of the conjugacy classes in a Sylow p-subgroup U(q) of a finite Chevalley group G(q), 
valid when q is a power of a good prime for G(q). In this paper we develop this algorithm 
and discuss an implementation in the computer algebra language GAP. Using the resulting 
computer program we are able to calculate the parametrization of the conjugacy classes 
in U(q), when G(q) is of rank at most 6. In these cases, we observe that the number of 
conjugacy classes of U (q) is given by a polynomial in q with integer coefficients. 



1. Introduction 

Let U„(g) be the subgroup of GL n (q) consisting of upper unitriangular matrices. A long- 
standing conjecture, attributed to G. Higman [12] states that the number of conjugacy classes 
of U n (g) is given by a polynomial in q with integer coefficients. This has been verified by 
computer calculation for n < 13 in the work of A. Vera-Lopez and J. M. Arregi, see [19]. 
This conjecture has generated a great deal of interest, see for example [17] and [18]. 

The equivalent problem of counting the number of (complex) irreducible characters of 
U n (g) has also attracted a lot of attention, see for example [15], [13] and [14]. Thanks to 
work of M. Isaacs [13], the degrees of the irreducible characters of U n (g) are known to all be 
powers q d of q and all exponents d occur for < d < fi(n), where fi(n) is an explicit upper 
bound due to work by G.I. Lehrer [15]. It was conjectured by Lehrer [15] that the number 
of irreducible characters of U n (g) of degree q d is a polynomial in q with integer coefficients 
only depending on n and d; this conjecture clearly implies Higman's conjecture. 

It is natural to consider the analogue of Higman's conjecture for other finite groups of Lie 
type. Below we introduce some notation in order to discuss this analogue for the case of 
finite Chevalley groups. 

Let G be a split simple algebraic group defined over the finite field W p of p elements, 
and assume that p is good for G. For a power q of p we write G(q) for the finite group 
of F g -rational points of G; this is a finite Chevalley group. Let U be a maximal unipotent 
subgroup of G defined over F p , so that U(q) is a Sylow p-subgroup of G(q). 

In this paper, we describe an algorithm that calculates a parametrization of the conjugacy 
classes of U(q). We have implemented this algorithm in the computer algebra language GAP 
[5]. The algorithm is based on the outline given by the first author in [8]. The output of the 
computer program allows one to calculate the number k(U(q)) of conjugacy classes of U(q). 
Using our computer program we have proved the following theorem. 
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Theorem 1.1. Let G be a split simple algebraic group defined over ¥ p , where p is good for 
G. Let U be a maximal unipotent subgroup of G defined over ¥ p . Let q be a power of p. 
If the rank of G is at most 6, then the number of conjugacy classes of U(q) is given by a 
polynomial in q with integer coefficients (and the polynomial itself is independent ofp). 

We have explicitly calculated the polynomial k(U(q)) for G(q) of rank at most 5; these 
polynomials are presented in Table 1 in Section 4. For G(q) of rank 6, the output of the 
computer program describing the parametrization of the conjugacy classes of U (q) is long 
and rather complicated. It is possible to check from this output in a short amount of time 
that k(U(q)) is given by a polynomial in q; however, a great deal of calculation would be 
required to explicitly compute this polynomial. We explain this in more detail in Section 4. 

Our computer program makes calculations using a Z-form of the Lie algebra u of U. In 
these calculations "implicit divisions" are made, which lead to the output not being valid 
for a finite number of primes p that are recorded within the program. In all case that we 
have run the program the only primes that are output are bad primes for G. A simple 
modification could insist on not dividing by certain prime numbers, and in this way it is 
possible to deal with the situation if any good primes were output. We note that the results 
on which our algorithm is based require the assumption that p is a good prime for G, so 
the output of our algorithm is not valid for bad primes even if there have been no implicit 
divisions. 

We have adapted the computer program so that it is also possible to use it to calculate 
the number of ?7(g)-conjugacy classes in certain subquotients of U(q). More precisely, let 
B = Ng{U), be the Borel subgroup of G corresponding to U, then for normal subgroups 
N D M of B that are contained in U, we can calculate a parametrization of the U(q)- 
conjugacy classes in the quotient N(q)/M(q). We have made a number of such calculations 
in case G is of rank greater than 6. In all cases where we have calculated the number of 
[/(g)-conjugacy classes in N(q)/M(q), we observe that it is given by a polynomial in q. In 
Table 2 in Section 4, we give the number k(U(q), U^ l \q)) of [/((^-conjugacy classes in the /th 
term U^ l \q) of the descending central series of U(q) for G of exceptional type and certain /. 

Generalizing a theorem of J. Alperin [1], the authors showed in [11, Thm. 4.6] that if the 
centre of G is connected, then the number k(U(q),G(q)) of conjugacy classes of U(q) in all 
of G(q) is a polynomial in q with integer coefficients (in case G(q) has a simple component 
of type E 8 , we require two polynomials depending on the congruence of q modulo 3). The 
theorem of Alperin [1] can be viewed as support for Higman's conjecture. Analogously, [11, 
Thm. 4.6] suggests that, for G not of type E$, the number of conjugacy classes of U(q) 
is given by a polynomial in q. The results of our computer calculations give supporting 
evidence for this behaviour, and we thus propose the following analogue and extension of 
Higman's conjecture for arbitrary finite Chevalley groups. 

Conjecture 1.2. Let G be a split simple algebraic group defined over ¥ p , where p is good 
for G. Let U be a maximal unipotent subgroup of G defined over ¥ p . Let q be a power of p. 
If G is not of type E 8 , then k(U(q)) is given by a polynomial in q with integer coefficients 
(and the polynomial itself is independent of p). If G is of type E 8 , then k(U(q)) is given by 
one of two polynomials depending on the congruence class of q modulo 3. 

The dependence of k(U(q)) on the congruence class of q modulo 3 in the E s case in 
Conjecture 1.2 is suggested by the E 8 case in [11, Thm. 4.5(h)]; though we do not wish to 
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rule out the possibility that there is just one polynomial. As indicated above, due to the 
complexity of the computation, it is not feasible to run our computer program in case G 
is of type E 8 . In fact, as shown in Table 2, at present we have only been able to calculate 
k(U(q), U^ l \q)) explicitly for / > 10; we have dim [J* 10 ) = 52 and dimf/ = 120 demonstrating 
the difficulty of running our program for E 8 . 

As mentioned above our algorithm is not valid for bad primes. It is possible to calculate 
the £/(g)-conjugacy classes for G of type B 2 and p = 2 by hand. In this case we have 
k(U(q)) = 5(q — l) 2 + 4(q — 1) + 1, which is a different polynomial to the one for good 
primes given in Table 1; this is due to degeneracies in the Chevalley commutator relations. 
In addition, [11, Thm. 4.6] is only valid for good primes, so we choose not to make any 
conjecture for bad primes. 

From our calculations we can observe that each polynomial k(U(q)), for G of rank 5 or less, 
when written as a polynomial in q— 1 has non-negative integer coefficients, see Table 1. For G 
of type A r and r < 12 this was already observed in the explicit results of Vera-Lopez-Arregi 
[19]. It would be interesting to have a geometric interpretation of this positivity behaviour. 
In Section 4, we give a reason why these positivity phenomena hold for the cases that we 
have calculated. This is done by analyzing the calculations made by the computer program. 
We expect that if Conjecture 1.2 is true, then the coefficients in k(U(q)) when written as a 
polynomial in q — 1 are always non-negative. 

In the cases where we have calculated k(U(q)), we have observed that k(U(q)) always has 
constant term equal to 1 when written as a polynomial in q — 1. In Section 4, we explain 
why this is necessarily the case whenever k(U(q)) is a polynomial in q. 

Another observation is that the polynomial k(U(q)) is the same for G of type B r and C r , 
for r = 3, 4, 5. It is likely that this is always the case for any r. We expect that this should 
be explained by the duality of the underlying root systems. 

Although k(U(q), G(q)) is always a polynomial in q (assuming that G has connected centre 
and taking into account that for G of type E 8 there are two polynomials) and our calculations 
here tell us that k(U(q)) is a polynomial for low rank, the number k(U(q), B(q)) of U(q)- 
conjugacy classes in B(q) is not always a polynomial in q. Indeed for G of type G 2 it is 
shown in [9, Exmp. 4.8] that the number k(B(q),U(q)) of 5(g)-conjugacy classes in U(q) 
is given by two polynomials: q + 15 if q is congruent to —1 mod 3; and q + 17 if q is 
congruent to 1 mod 3 (we are assuming p is good for G). A general argument considering 
the commuting variety C(B,U) = {(b,u) G B x U \ bu = ub} shows that we always have 
k(U(q), B{q)) = {q- l) 2 k{B{q), U(q)). Thus k(U(q), B{q)) is not a polynomial in q for G of 
type G 2 . 

Our algorithm calculates a family of varieties X c that parameterize the conjugacy classes of 
U ; moreover, these varieties are defined over F p . The algorithm determines the polynomials 

defining the X c as locally closed subsets of (F*) mc for certain m c e Z> . The varieties X c 
are determined with a backtrack algorithm using a depth- first search. The conjugacy classes 
of U (q) are parameterized by the F ? -rational points of the varieties X c and it is possible to 
count these points. 

The idea behind the algorithm is similar to that for the algorithm used by Biirgstein and 
Hesselink in [2] for calculating the adjoint orbits of B in u = Lie U ; we remark that the 
algorithm in loc. cit. was not written to give a complete description of the £?-orbits in u. 
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In addition, our algorithm generalizes the one used in the work of Vera-Lopez and Arregi 
for the type A situation, see for example [19]. Finally, we remark that the algorithm of this 
paper uses ideas from the computer program described in [6] in previous work of the first 
author. 

We now give a brief outline of the structure of this paper. In Section 2, we introduce the 
notation that we require and recall the relevant results from [8] and [9]. Then in Section 3 
we describe the algorithm and its implementation in GAP. Finally, in Section 4 we discuss 
the results of our computations. In particular, we present explicit values for k(U(q)) for G 
of rank at most 5 (Table 1) and the values of k(U(q), U®(q)) for some cases where G is of 
exceptional type (Table 2). 

As general references for algebraic groups defined over finite fields we refer the reader to 
the books by Carter [3] and Digne-Michel [4]. 

2. Notation and recollection 

Let p be a prime and let G be a split simple algebraic group defined over the finite field of 
p elements F p . We assume throughout that p is good for G and we write k for the algebraic 
closure of F p . 

Fix a split maximal torus T of G and let $ be the root system of G with respect to T. For 
a root a G $ we choose a generator e a for the corresponding root subspace Q a of q = LieG. 
Let B D T be a Borel subgroup of G that is defined over F p . Let U be the unipotent radical 
of B and let u = Lie U. Let <3> + be the system of positive roots of $ determined by B. The 
partial order on <3> determined by $ + is denoted by -<. 

For a power q of p, we write G(q) and U (q) for the finite groups of Fq-rational points of 
G and U respectively. We write u(q) for the Lie algebra of F g -rational points of u. 

We now recall some results from [8] and [9] on which our algorithm for calculating the 
conjugacy classes of U(q) is based. Thanks to [7, Thm. 1.1], there are generalizations of 
some results in [8], as explained in [7, §6]; below we state the general versions without 
further comment. For the remainder of this section we fix a power q of p. 

The theory of Springer isomorphisms can be used to show that the conjugacy classes of 
U(q) are in bijective correspondence with the adjoint [/(g)-orbits in u(q), see for example 
[8, Prop. 6.2]. For the remainder of this paper, we will consider the adjoint orbits of U in u 
rather than the conjugacy classes of U, as this is more convenient for our purposes. 

Next we discuss the notion of minimal representatives of [/-orbits in u, and how they are 
used to partition the set of [/-orbits in u. The reader is referred to [8, §5 and §6] and [9, §3 
and §4] for full details. 

We fix an enumeration of the set of positive roots <£> + = . . . ,@n}, such that % < j 
whenever ^ f3j, and define the sequence of B-sub modules 

u = m 2 ■ ■ • 2 mjv = {0} 

of u by rrij = J2f=i+iS/3j- We consider the action of U on successive quotients Uj = u/rtij 
induced from the adjoint action of U on u. We note that the parametrization of the adjoint 
[/-orbits described below depends on the choice of the enumeration of $ + . 
Let x G u and consider the set 

x + ke^ + ttij = {x + Ae ft + rtij | A G k} C u». 
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By [8, Lem. 5.1], for x G u either: 

(I) all elements of x + ke^ + are [/-conjugate; or 
(R) no two elements of x + ke^ + rrij are [/-conjugate. 

We say that 

• i is an inert point of x if (I) holds; 

• i is a ramification point of x if (R) holds. 

An element x + trij = X^=i x j e Pj + m « °f u « ls sa ^ to the minimal representative of its 
[/-orbit in Uj if = whenever j is an inert point of x. It follows from [8, Prop. 5.4 and 
Lem. 5.5] that each [/-orbit in u« contains a unique minimal representative; in particular, 
this holds for the action of U on u. 

Thanks to [8, Prop. 4.2 and Lem. 5.7], we have that % is an inert point of x G u if and 
only if dim c u (x + ttij) = dim c u (:r + rrij_i) — 1; if i is a ramification point of x, then we have 
dim c u (x + rrij) = dim c u (x + ntj_i). Here c u (x + rrij) is the centralizer of x + ttij for the action 
of u on Uj induced from the adjoint action of u on itself. 

The above discussion implies that the adjoint orbits of U in u are parameterized by their 
minimal representatives. Further, the set of minimal representatives can be partitioned into 
sets X c for c G {I, R}^: the set X c is defined to consist of the minimal representatives x G u 
of the [/-orbits in u such that for alH = 1, . . . , N we have that % is an inert point of x if and 
only if q = I. Thanks to [9, Prop. 2.4], each of the sets X c is a locally closed subset of u, 
and therefore has the structure of an algebraic variety. 

The above partition of the [/-orbits in u can be refined to be indexed by iV-tuples c G 
{I, R n , Ro}^ as follows. For c G {I, R n , Ro}^, the set X c is defined to consist of the minimal 
representatives x = Yl x i e i3i G u of the [/-orbits in u such that for all % = 1,...,N we 
have that i is an inert point of x if and only if q = I; and if q ^ I, then Xi — if and 
only if Ci = R . Thanks to [9, Lem. 4.2], each of the sets X c is a locally closed subset of 
u, and therefore has the structure of an algebraic variety. In fact, X c is a subvariety of 
{(^■) c . =Rn | x 3 G k*} = {k x ) m % where m c = \{j | Cj = R n }|- 

We now explain how the above parametrization of the [/-orbits in u descends to give a 
parametrization of the U (g)-orbits in u(q). The reader is referred to [8, §6] for further details. 

Thanks to [8, Prop. 4.5], we have that for x G u the centralizer Cjj(x) of x in U is connected. 
This implies that the [/(g)-orbits in u(q) correspond bijectively to the [/-orbits in u that are 
defined over ¥ q . Let x G u be the minimal representative of its [/-orbit. Then, by [8, Lem. 
6.3], the orbit U ■ x is defined over ¥ q if and only if x G u(q). We require that the definition 
of G over W q is split for this last assertion, and this is the reason for this assumption. If the 
definition of G over ¥ q is not split, then it is a non-trivial task to determine which minimal 
representatives correspond to orbits defined over ¥ q . 

It follows from the above discussion that the adjoint orbits of U(q) in u(q) are parameter- 
ized by the minimal representatives of the [/-orbits in u that lie in u(q). In turn these minimal 
representatives are partitioned into the sets X c (q) of F^-rational points of the varieties X c , 
for cG {I, R n ,R } w . 

3. The algorithm 

In this section we develop the algorithm outlined in [8] for calculating the parametrization 
of the adjoint [/-orbits in u. The idea is to calculate the polynomials defining the varieties X c 
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for c G {I, R n , Ro}^ as locally closed subsets of (k x ) mc . We present the algorithm, and then 
explain why the algorithm does indeed calculate a parameterization of the adjoint [/-orbits 
in u. Next we discuss two modifications that are used in the GAP implementation, before 
briefly explaining the implementation. Finally, we explain how the output of the computer 
program is used to calculate k(U(q)). 

In order to explain the algorithm we have to introduce some more notation; we continue 
to use the notation given in the previous section. 

We wish to consider all primes p simultaneously, so we need a Z-form of $j. Let Qc be the 
complex simple Lie algebra of the same type as g. Fix a Chevalley basis of 0c and let Qi be 
the corresponding Z-form of g c . We let 



where Z[ti, . . . ,t m ] is the polynomial ring in m indeterminates t±, . . . ,t m . We denote by 
ep 1 , . . . , Cj3 N the elements of the Chevalley basis of Qz corresponding to $ + , which is enumer- 
ated as in the previous section. These elements form a Z-basis of uz, and by a minor abuse 
we also consider them as elements of both u and u. 

Let c G {I,R n ,R } 1 for some % < N. For j = 1, . . . ,m c , we define (3 c j G $ + by setting 
Pc,j = Pi, where / is the jth smallest element in {m \ c m = R n }. We associate to c the 
element 



j'=i 

Given r = (r 1? . . . , r m J G k mc , we write x c {r) for the element of u obtained by substituting 
tj = Tj in x c (t), i.e. 



The variety X c is a locally closed subset of {x c (t) \ r G (k x ) mc } = (k x ) mc . Therefore, 
there are subsets A l c) . . . , A l £ and B\, . . . , B l c c of k[t 1: . . . , t m J C k[t 1: . . . , t m ] such that X c is 
the disjoint union of the sets 



for % — 1, . . . , l c . In fact the polynomials in the sets A l c and B l c can be taken to have integer 
coefficients; this is due to the integrality of the Chevalley commutator relations. The purpose 
of our algorithm is to determine certain choices for the sets A % c and B l c . 

We note here that it is most often the case that we can take l c = 1 and A\ = B\ = 0. 
The values of c for which this is not the case in some sense explain why the determination 
of the conjugacy classes of U(q) is complicated in general. We also remark that it is often 
the case that X c = 0, which corresponds to the case / = 0. 

We now introduce some notation needed in order to say how the sets A l c and B l c are 
determined in the algorithm. Let yi, . . . ,Un be variables. We may write: 



m := max{m c | c G {I, R n , Ro} , X c ^ 0}, 
where m c = \{i \ Ci = R n }|, as defined earlier. We define 

u = u z <8> z Z[ti,...,t m ], 





X l c = {x c (t) I /(r) = for all / G A\ and g(r) ^ for all g G 



f N INN 



U=l J 3=1 k=l 
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where each Pj k (t) G Z[£i, . . . , t m ] is linear: this is easily achieved using the Chevalley com- 
mutator relations for Qc- R is then the case that dimc u (a; c (r) + ttij) is the dimension of the 
solution space of the system of linear equations: 

N 

k=i 

for j = l,...,i. 

We are now in a position to describe our algorithm. It calculates sets of polynomials 
that determine certain choices of the varieties X l c . These polynomials are calculated using 
the Z-form Ui of u and they have rational coefficients. Observe that in step (2) (a) of the 
algorithm there may be "implicit divisions" by certain primes for which the output will not 
be valid. This is discussed in more detail later in this section. 

The algorithm uses a backtrack algorithm with a depth-first search to calculate certain 
choices of sets A % c and B l c that determine the varieties X l c for each c G {I, R n , Ro}^, and 
% — 1, . . . , l c . In the algorithm we require a total order on Z[£i, . . . , t m ]; we use the order 
defined in precedence by the number of terms, total degree and leading coefficient (with 
respect to the degree then lexicographic order on monomials). 

The variables used in the algorithm are: 

• the "current string" c is an element of {I, R n , R }* for some % and determines x c (t) G u 
as above; 

• the set of "satisfied" polynomials A is a subset of Z[ti, . . . ,t m ]; 

• the set of "unsatisfied" polynomials B is a subset of Z[ii, . . . ,t m ]; 

• the matrix Q(t) is an element of Matj j jv(Z[ii, . . . ,t m ]), which is obtained from the 
matrix (Pjf.(t)) by "row reducing" the first % rows; 

• the "pivot string" ir is an element of {0, 1, . . . , N} 1 , which "records the columns used 
in the row reductions" ; 

• the stack S = {(c, A, B, 7r, Q(t))} is an (ordered) subset of 

N N 

(J{I, R n , Ro} 1 x P(Z[ti, . . . , t m ]f x |J{0, 1, . . . , NY x Matjv(Z[ti, . . . , t m }), 

i=l i=l 

which contains variables to be considered later in the algorithm; and 

• the output set O is a subset of {I, R n , Ro}^ x P(Z[ti, . . . , t m ]) 2 . 

Here, P denotes the power set. 

The stack is required to be ordered as the algorithm takes elements from the "top" of the 
stack. The element at the top is the one that has been most recently added and is denoted 
by top (S). The initial configuration in the algorithm is as follows: 

• c := (R n ); 

• A := 0; 

• B := 0; 

• vr := (0); 

• Q(t) is the 1 x N matrix with all entries equal to 0; 

• S:= {(R , 0,0,(0), Q(t))Y and 

• O := 0. 
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Now we explain the next step in the algorithm; we have numbered the steps in the algo- 
rithm, so that we can refer back to it in the explanation given afterwards. 

(1) If the length of c is N, then we are finished with this string. We set: 

(a) O :=OU{(c,A,B)} 
If S — 0, then we finish. 

Else we make the following changes to the variables: 

(b) (c, A, B, Ti, Q) := top(S'); and 

(c) S:=S\{top(S)}. 

(2) If the length of c is i — 1 < N, then we proceed by making the iih row reduction for 
the matrix (Pj k (t)) as defined above. Note that Q(t) is the matrix resulting from the 
first i — 1 row reductions. We first append the ith row of (P? fc (£)) to Q(t) and then 
make the row reduction as follows: 

(a) for j = 1, . . . , i - 1, if TTj ^ we set Q^t) := Q^. (*)&(*) - QJ i)r .(f )<?;(*), where 
Q'j % .{t) is Qj,TTj(t) divided by the highest common factor of Qj tWj (t) and Qi, nj (t), 
and Q' ij1T (t) is defined analogously. 

Let Ri be the set of non-zero polynomials in the ith row of Q(t) that are not divisible 
by any element of A. We next consider three cases: 

(b) Ri = 0. We update the variables as follows: 

(i) 7T := (tt, 0) 

(ii) c := (c, R n ); and 

(iii) S:=SU{((c,R ),A,B,7r,Q)}. 

(c) Ri ^ and there is some non-zero element of Ri that is a monomial or divides 
some element of B. We let Qi t i(t) be the least such polynomial with respect 
to our chosen order on the set of all polynomials. We update the variables as 
follows: 

(iv) c := (c, I); and 

(v) tt := (tt,/). 

(d) Otherwise, we pick a least element Qi,i(t) of Ri. We update the variables as 
follows: 

(vi) S:=SU{(c,AU{Q hl (t)},B,n,Q)}; 

(vii) c := (c, I); 

(viii) B := BU{Qi,,(t)}; and 
(ix) tt := (tt,/). 

The output of the algorithm is a collection of triples (c,A,B). Each of these triples 
determines a subvariety 

(3.1) X M)B =J J] r^ Ci . G X c | /(r) = for aJl / G A and ^(r) ^ for all ^ G S 

^ j-.Cj=R n 

of X c . For fixed c, we can write A\, . . . , A[, c and B^, . . . , B l c c for the sets A and i? occurring 
in a triple (c, A, B). These are the determined choice of the sets A l c and B l c such that X c is 
the disjoint union of the varieties X\ as defined earlier; each X l c is equal to the corresponding 

We now explain why the output of our algorithm can be used to determine all the minimal 
representatives of the £7(g)-orbits in u(q), for almost all primes p. As explained above there 
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may be "implicit divisions" in step (2) (a) of the algorithm leading to a finite number of primes 
for which we cannot determine the minimal representatives of the [/(g)-orbits in u(q). We 
choose not to give a formal proof of the correctness of the algorithm as this would be very 
technical and just give an outline. We argue by induction on % to show that the algorithm 
determines varieties X c a,b,i (as defined below) from which one can calculate all the minimal 
representatives of the [/(^-orbits in Ui(q) for each i (and valid p). We do not discuss the 
part of the algorithm dealing with the row reductions in (2) (a), as this is elementary. 

In the discussion below we often speak of all relevant r G k m . When considering the triple 
(c, A, B), this means all r G k m such that /(r) = for all f e A and g(r) ^ for all g G B, 
i.e. all t G k rn for which x c (t) G X Cj a,b- 

We need to explain our inductive hypothesis as this is not compatible with the depth- 
first search used in the algorithm. To do this we must consider all triples (c, A, B) with 
c of length i that occur during the running of the algorithm. For such (c, A, B) varieties, 
X c ,a,b can be defined as in (3.1). Then the inductive hypothesis says that the varieties 
X c ,A,B,i = {x + trij | x G X c a,b} gi ye all minimal representatives of the [/-orbits in u«. 

So assume the inductive hypothesis for i — 1. Then we have to show that for each (c, A, B) 
with c of length i and all minimal representatives of [/-orbits in Uj of the form x c {r) +Ae / g i +rrij 
with x c (t) G X c> a,b and A G k, we have that x c {r) + Ae^ lies in some variety of the form 
X c ' t A',B' where d = (c, Z) and Z G {I, Ro, R n }- To do this we have to consider the steps in 
(2) in the algorithm. 

After the row reduction made in (2) (a), we have the row reduced matrix Q(t) and the 
set Ri. As explained earlier we have that i is an inert point of x c {r) G X c a,b if and only 
if dimc u (a; c (r) + ntj) = dimc u (x c (r) + rrij_i) — 1. Also the dimension of c u (x c (t) + rtij) is 
the rank of the row reduced matrix Q{r). Therefore, % is an inert point of x c (r) if and only 
if f(r) 7^ for some / G Ri. Note that we only have to consider the polynomials in R4, 
because for any non-zero entry in the ith row of Q(t) that is divisible by some polynomial 
f(y) G B, we automatically have /(r) = for all relevant r. At this stage one ideally wants 
to determine for which values of r there is some / G R4 with /(r) 7^ 0. However, this is a 
difficult task if Ri contains several polynomials. Also we require a fixed value for the next 
entry of 7r, so that we can perform the row reductions later in the algorithm. So we proceed 
by considering the three cases (2)(b)-(d) in the algorithm: 

(b) In this case it is clear that % is a ramification point of x c (r) for all relevant r. There- 
fore, x c (t) + Xefc + rrij is the minimal representative of its [/-orbit in Uj for all A G k. 
We see that the strings d corresponding to all such minimal representatives are 
passed on in the program: the non-zero values of A correspond to the updated strings 
d = (c, R n ) in (ii), and the case A = corresponds to the string d = (c, Ro) added to 
the stack. 

(c) For this case, it is clear that f(r) 7^ for all relevant r, where f(t) = Qi,i(t) is the 
least element of Ri that is either a monomial or divides some element of B. Therefore, 
% is an inert point of x c (r) for all relevant r. Thus the only minimal representative 
of its [/-orbit in u, of the form x c (t) + Ae^. + rrij is when A = 0. This minimal 
representative corresponds to the updated string d = (c, I) in (iv). 

(d) This case is more complicated. We consider the least element f(t) = Qi : i(t) of Ri 
and the following cases for r: 
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(I) If /(t) 7^ 0, then % is an inert point for x c (r). For such r, the only minimal 
representative of its [/-orbit in Uj of the form x c (t) + Ae^ + rrij is when A = 0. 
This minimal representative corresponds to the updated string d = (c, I) in (vii), 
along with / being added to B in (viii). 
(II) If /(r) = 0, then we cannot say whether % is an inert or ramification point of 
x c (t). The element (c, A U {/}, B, ir, Q(t)) added to the stack in (vi) will be 
considered later in the algorithm. For this case Ri will have fewer elements, 
and will be considered again either in case (a) or (c). This process will finish, 
and it is straightforward to see that when this happens all the required triples 
(d, A', B') with d = (c, Z), will have occurred. 

Putting together the above case analysis, verifies the inductive step. Therefore, for all 
minimal representatives x c (r) of the [/-orbits in u, a corresponding triple (c, A, B) occurs at 
some point of the program. This triple is added to the output in (l)(a), which means that 
the output of the algorithm determines all minimal representatives, as claimed. 

Next we discuss two modifications to the algorithm that we make for its implementation 
in GAP. These changes are made in order to speed up the computations; we chose not to 
include them in the above description of the algorithm for simplicity. 

Our first modification allows us to reduce the complexity by using the action of the maxi- 
mal torus T to "normalize" certain coefficients to be equal to 1. Let c G {I, R n , Ro}^ and let 
x c{t) — 'YlT=i^i e Pci be defined as above. Suppose that {(3 c ,i | 2 G J} is linearly independent 
for some J C {1, . . . , m c }. Then for every r G (k x ) mc there is some a G (k x ) mc with Oi — 1 
for alH G J and x c (r) is conjugate to x c (a) via T, i.e. the action of T can be used to "normal- 
ize Tj = 1" for all % G J. In the computer program we replace x c (t) by Yli^j e Pc % + ^2i<£j ^ e /? c * 
thereby reducing the number of indeterminates required. This speeds up the row reduction 
of the matrix Q(t) significantly. Some care is needed when using this modification, as the 
centralizer C t (x c (t)) can be disconnected. If this is the case then it can become difficult to 
determine the elements of X c (q) from the "normalized" elements of X c . This problem can 
be resolved by not allowing certain normalizations; we omit the technical details here. 

The second adaptation deals with "easy" elements of the set A. If there is a linear 
polynomial in A, then we may simplify future checks by "making a substitution" . If ti—a(t) G 
A, where i G {1, . . . , m c } and a(t) is a linear polynomial not involving ti, then we can make 
the substitution ti = a(t) in x c (t), in the polynomials in A and B, and in the matrix Q(t); 
we then remove ti — a(t) from A. This modification reduces the number of indeterminates, 
and also the number of polynomials in A. This helps to speed up the program. 

We next explain a check that has to be included in the program to see which primes 
the output is valid for. When making the row reductions in (2) of our algorithm, there 
are implicit divisions by certain integers. Essentially, we need to be able to divide by the 
polynomials in the set B in (2) (a). Therefore, we keep track of the primes dividing their 
leading coefficients, and these primes are output by the program. The output of our program 
cannot be used to determine the minimal representatives of the U (g)-orbits in u(q) for such p. 
The only primes output by the program for the cases that we have calculated are bad primes 
for G. If a good prime p were output by the program, then it would be straightforward 
to adapt the program to insist that there are no implicit divisions by p, and running this 
modification would give an output valid for p. 
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The algorithm is implemented in GAP with the two modifications and the check for primes. 
This is achieved using the functions for Lie algebras and polynomial rings in GAP. This allows 
us to define u within GAP and therefore allows us to obtain the matrices (Pj k (t)) that we 
row reduce using the method given in (2) (a) of the algorithm. The implementation is based 
on the algorithm given and the two modifications discussed above. We choose not to include 
any of the technical details. 

In the next section we present the values of k(U(q)) that we have calculated from the 
output of our program. Each of the varieties X c ^ a ,b is defined by polynomials with integer 
coefficients, so is defined over ¥ p . We have 

Hu(q))= J2 \ x caAq)\- 

(c,A,B) 

We can, therefore, calculate k{U(q)) by calculating \X Cj a,b{c[)\ for all triples (c, A,B). If the 
polynomials in A and B are not too complicated, then this can be achieved quite easily. We 
discuss this below. 

It is most commonly the case that c occurs in just one triple (c, A, B) for which both A 
and B are empty. In which case it is easily seen that X c = X C) a,b and |X c (g)| = (q — l) mc . 
The next simplest case is when AU B has one element that is linear. For example, consider 
the polynomial t x - 1: if A = {ti - 1} and B — 0, then \X CtAjB (q)\ = (q - l) mc_1 ; and if 
A = and B = {t ± - 1}, then \X c>AB (q)\ = (q - l)"^ 1 ^ - 2). More complicated sets A 
and B that we need to consider require a little thought to calculate \X CtAiB (q)\. 

As the rank of G increases the polynomials become more complicated. For the F4, B5 and 
C5 cases we get a number of quadratic polynomials. For the rank 6 cases, the polynomials 
become more complicated still and the number of triples (c, A, B) with A or B non-empty 
gets large. From the output of the program we can view all the polynomials occurring in the 
sets A and B. One can check that for each of these polynomials it is possible to "solve for 
one indeterminate in terms of the others". With a little further consideration one can see 
that this means each of the sets X CjAB (q) has size a polynomial in q. However, the number 
of such X c> a,b(q) is so large that it would be rather time consuming to calculate k(U(q)) 
explicitly. 

4. Results 

In this final section we present some explicit results of our computations and go on to 
discuss some interesting features of the output. 

In Table 1 below we present the polynomials k(U(q)) for G(q) of rank at most 5; in this 
table we let v — q — 1 to save space. We include the values for G of type A r for completeness, 
though these polynomials have been known for some time, thanks to the work of Vera-Lopez 
and Arregi referred to in the introduction. Also, as discussed below, the value of k(U(q)) is 
the same for G of type B r and C r , so we only include this polynomial once. 

We make some comments about the polynomials in Table 1. We start by making the 
observation that k(U(q)) considered as a polynomial in v — q—1 has non-negative coefficients. 
For the case G is of type A r (r < 12), this was observed by Vera-Lopez and Arregi in [19]. 
It would be interesting to have a geometric explanation of these positivity phenomena. 

We give a heuristic idea why this occurs for the cases that we have calculated by considering 
the partition of the conjugacy classes used in our algorithm. As discussed at the end of the 
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G 


k(U(q)) 




v 4- 1 




7 ;2 _l_ Q71 -I- 1 


Bo 






^3 _|_ 5^2 _|_ _|_ 1 


A? 


2v 3 4- 7t> 2 4- 6v 4- 1 

til/ | 1 \J | v_/ 1/ | -L 




7 ,4 4. fi ?; 3 1 lfi,,2 1 q ? . 1 1 

(y | uu | _l \j u | kJ \J \ J- 




O U \^ ZjW V \^ ZiO U \^ 1U U \^ X 


BaICa 


1; 6 1 1 1 7 .5 1 4 o ? .4 1 oo ? .3 1 C4 ?; 2 , 1 c , 1 


Da 


2-j; 5 4- 157, 4 + 36t; 3 4- 34t; 2 4- \2v 4- 1 


Fa 
A 4 


7,8 i 0^7 — 1 — 4q ? ,6 — 1 — i24 v 5 + 25Qv 4 4- 288v 3 4- 140v 2 4- 24?; 4- 1 


A 5 


4- 18w 5 4- 70v 4 + I05v 3 + 65v 2 + 15v + 1 


-B5/C5 


2w 8 4- 24w 7 4- 132tj 6 4- 395v 5 + 630v 4 + 500v 3 + 180v 2 4-25^ + 1 


D 5 


2v 7 + 22v 6 + 106v 5 + 235v 4 + 240v 3 + 110v 2 + 20v + 1 



TABLE 1. k(U(q)), as polynomials in v — q — 1 



previous section, the number of the F g -rational points of the varieties X l c is most commonly 
|X*(g)| = v mc . Although there are some values of c and % for which |X*(g)| is a polynomial 
in v with negative coefficients, these negative coefficients are few enough so that they are 
cancelled by the families of size v mc . 

We observe that the constant coefficient in k(U (q)) as a polynomial in v is always 1. This is 
explained by the action of the split maximal torus T of G on each X c for all c G {I, R n , Ro}^- 
This action is non-trivial unless q = Ro for all i, so that X c = {0}. It is easy to see that if 
X c ^ {0}, then the orbits of T(q) on X c (q) are all of size v a /b for some a,b e Z>i, so |X c (g)| 
is divisible by t> a /6. This implies that the constant coefficient in k(U(q)) as a polynomial in 
v must be 1 (corresponding to the zero orbit). 

We now comment on the fact that the value of k(U(q)) is the same for G of type B r and 
C r , for r = 3, 4, 5. One can see that the groups U (q) are not isomorphic for G of type B r and 
C r : thanks to a result of A. Mal'cev [16], the maximal size of an abelian subgroup of U(q) is 
different for G of type B r and C r . Using the variation of our program discussed below, one 
can also show that the number of £/(g)-conjugacy classes in the derived subgroup U^^q) 
of U(q) are different for G of types B r and C r , for r = 3,4,5. It would be interesting to 
have a reason for the coincidences in the numbers k(U(q)); we expect it should be explained 
by the duality of the root systems of type B r and C r , see for example [3, Ch. 4] for similar 
phenomena. 

12 



As mentioned in the introduction, we have adapted our program to consider the action of 
U on certain subquotients M/N. The adaption is valid when M D N are normal subgroups 
of B contained in U. The algorithm runs in essentially the same way: one has to replace the 
filtration of u by an analogous filtration of m/n, then change the initial configuration and 
the point at which variables are added to the output set O accordingly. 

In Table 2 below we give some values of k(U(q),U^(q)) for G of exceptional type. We 
recall that the descending central series of U is defined by = U and = \U^ l ~ 1 \U\ 
for / > 1. The cases that we have included are those for which we are able to compute 
k(U(q), U®(q)) in a reasonable amount of time and for which there is an infinite number of 
.B-orbits in = LieC/®; we refer the reader to [10] for a classification of all cases when 
there is only a finite number of 5-orbits in for G of exceptional type. 



G 


I 


k(U(q),U«\q)) 


F 4 


1 


v 7 + 7v 6 + 2Av 5 + Q3v A + 119v 3 + 88v 2 + 20v + 1 




2 


2v 5 + Uv 4 + 50v 3 + 58v 2 + 17^ + 1 




3 


2v A + 18v 3 + 35v 2 + 14^ + 1 


Eq 


1 


v w + I0v 9 + A7v 8 + I53v 7 + 435w 6 + 993v 5 + 1315v 4 + 868v 3 + 255v 2 + 30^ + 1 




2 


2v 7 + 28v & + 160w 5 + 38Qv A + 404w 3 + IQhv 2 + 25^ + 1 




3 


v 6 + llv 5 + 70v 4 + U8v 3 + 95v 2 + 20v + 1 


E 7 


4 


v 9 + 13v 8 + 9Av 7 + 512v 6 + WOOv 5 + 2312w 4 + U99v 3 + 395v 2 + 38v + 1 




5 


v 8 + 10v 7 + Q3v 6 + 292v 5 + Q85v 4 + 700v 3 + 2Q0v 2 + 32^ + 1 




6 


3v & + 39v b + \72v A + 3l2v 3 + 170v 2 + 27^ + 1 


E 8 


10 


v 9 + 17v 8 + 135v 7 + 7l9v 6 + 2568w 5 + 4652w 4 + 3014w 3 + 699v 2 + 52^ + 1 




11 


v 8 + 12v 7 + 92v 6 + 518v 5 + 1766v 4 + 1693w 3 + 516v 2 + 46v + 1 




12 


5v 6 + 67v 5 + 660v 4 + 96Av 3 + 386v 2 + 41u + 1 



Table 2. k(U(q), U {l \q)), as polynomials in v — q - 1 
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